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A multi-dimensional semiclassical description of excita- 
tion of a Rydberg electron by half-cycle pulses is developed 
and applied to the study of energy- and angle-resolved ioniza- 
tion spectra. Characteristic novel phenomena observable in 
these spectra such as interference oscillations and semiclas- 
sical glory and rainbow scattering are discussed and related 
to the underlying classical dynamics of the Rydberg electron. 
Modifications to the predictions of the impulse approximation 
are examined that arise due to finite pulse durations. 

PACS numbers: 32.80.Rm, 03.65.Sq 

High power, nearly unipolar electromagnetic field 
pulses are a useful new spectroscopic tool in particu- 
lar for studying the dynamics of weakly bound Ryd- 
berg electrons. Various experimental and theoret- 
ical m investigations on the interaction between Ryd- 
berg atoms and HCPs have already been performed. In 
particular, the theoretical approaches which have been 
used so far involve cither fully quantum mechanical cal- 
culations, purely classical simulations, or simplified one- 
dimensional model problems. Fully quantum mechanical 
treatments are capable of yielding quantitatively exact 
results but due to the large spatial extension of highly ex- 
cited Rydberg states they demand a huge computational 
effort. Furthermore, they allow only for a restricted qual- 
itative understanding of the system, in general. Purely 
classical approaches are of interest as they yield insight 
into the underlying classical aspects of the motion of the 
Rydberg electron. But they cannot deal properly with 
quantum mechanical interference effects and are there- 
fore also of limited applicability. One-dimensional model 
problems, finally, are not capable of describing all spa- 
tial aspects of the excitation process. Thus, the natu- 
ral question arises whether a proper three-dimensional 
theoretical description of the excitation process can be 
developed which is numerically highly accurate even for 
very large principal quantum numbers, which gives di- 
rect insight into the underlying classical aspects of the 
dynamics of the Rydberg electron and which is capable 
of dealing with all quantum mechanical interference phe- 
nomena properly. 

In this Communication it is shown that such a the- 
oretical description can be developed on the basis of a 
multidimensional semiclassical description of the excita- 
tion process. This approach is easily applicable to any 



spatial or temporal pulse form of the exciting HCP and 
it is particularly accurate in the region of high princi- 
pal quantum numbers which is difficult to access by fully 
quantum mechanical calculations. Based on this treat- 
ment novel oscillatory structures are discussed which ap- 
pear in energy- and angle-resolved ionization spectra and 
which should be accessible to experimental observation 
in view of the recently developed imaging method || . It 
is shown that these structures are caused by quantum 
mechanical interferences between probability amplitudes 
that can be associated with classical trajectories of the 
ionized Rydberg electron. In order to demonstrate the 
numerical accuracy of the multidimensional semiclassi- 
cal approach a comparison with numerical results is pre- 
sented in the impulse approximation. Modifications of 
these oscillatory structures due to finite pulse durations 
are also discussed. 

Let us consider a typical HCP-excitation process. Ini- 
tially an atom is prepared in a Rydberg state |noZo TO o) 
with a large principal quantum number % > 1 and 
a small angular momentum quantum number <C n o- 
Around t = a linearly polarized HCP with vector po- 
tential ^4(x, t)e z and duration r interacts with the Ryd- 
berg electron. The time evolution of the state ^(x, t) of 
the Rydberg electron is determined by the Schrodinger 
equation 

#(x, t) = {i[- zV x - A(x, t)e z } 2 - Vfc(x)}V(x, t) (1) 

with the effective electronic potential Vc(x) (Hartree 
atomic units with e = fi = m c = 1 are used). As- 
suming axial symmetry of A(x, t) around the z-axis the 
z-component of the electronic angular momentum L z is 
conserved. As long as no ^ 1 the solution of Eq. (|l|) can 
be obtained to a good degree of approximation with the 
help of semiclassical methods. Thereby, one starts from 
the semiclassical expression for the initial Rydberg state 

(xoMomo) = A cl (x ){e^ r °^-*^ + c.c.}. (2) 
The amplitude of this state is given by Q] 

y, mo (e o ,0o) 

^d(xo) = ~ , = (3) 

r (n - a) 3 / 2 ^/27rp(r , e ) 

with the spherical harmonic Y™° , the radial momen- 
tum p(ro,eo) — -\/2(eo + l/ r o)i an d with the initial en- 
ergy eo = — [2(no — a) 2 ] -1 {a denotes the quantum 
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defect of the Rydberg electron). The phases of this 
state are determined by the classical eikonal iSo(?*o, e o) = 
J r dr'p(r', eo) — (Zo + 1/2)tt + ira. Equation (||) is valid 
for distances ro — |xo| of the electron from the nucleus 
which arc located outside the core region and well inside 
the classically allowed region. Semiclassically one has to 
associate with this particular initial state a bi- valued field 
of classical momenta p^^Xo) = ±V Xo 5'o(r , e ) which 
defines a Lagrangian manifold. Generalizations to other 
types of initial states which lead to other Lagrangian 
manifolds are possible according to the general frame- 
work of multidimensional semiclassical approximations 
||. Thus under the influence of the HCP two classical 
trajectories x±(i;xo) emanate from each point xo. They 
are solutions of the classical equations of motion of the 
Hamilton function H = [p — A(x, t)] 2 /2 + Vc(x) with ini- 
tial conditions [xo, p^^xo)]. From Eq. ([j]) and Maslov's 
generalization of the Van-Vleck propagator Q one ob- 
tains the probability amplitude of observing the ionized 
Rydberg electron with final asymptotic momentum p^, 
namely 

{^m^oo)) = Y,pf ) ^ [S]+w ^ )) ^ H ' 2] - (4) 



As this result is understandable in a straightforward way 
on physical terms the details of its derivation will be pre- 
sented elsewhere Q. 

The ionization amplitude of Eq. (|4|) is expressed 
as a sum of contributions of all classical trajectories 
j of the Rydberg electron with initial positions xj, 
which yield the asymptotic momentum p^ = p±(t — > 
oo; xo, p^^Xo)). The contribution of trajectory j to the 
classical angle- and energy-resolved ionization amplitude 
is determined by 



3 



dp [ J ] A d P { y f) A dpi n 



dxo A dyo A dzo 



-1/2 



(•5) 



The phases associated with these trajectories are deter- 
mined by the Morse indices Hj , by the classical actions 
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dt Xj (t) 



dp 



U) 



dt 



v (j) . n U) 



(6) 



and by the classical actions of the initial state Wj(x ~'' ) ) = 



±[S , (r j) , e ) - tt/4]. For Wj(-x.^') one has to choose the 
plus or the minus sign depending on whether the initial 
radial momentum of the classical trajectory j is positive 
or negative. The Morse index is equal to the number 



of sign changes of — ^ — . . , = Det — w , * — 

b ° dx f\dy a f\dza d(xa,ya,z a ) 

along the trajectory times their multiplicities H. In 
terms of the probability amplitude of Eq. (|J) the semi- 
classical angle- and energy resolved ionization probability 
is given by 



d 3 P 

u 1 ion 

de A dfl 



(pW\tp(t-> oo)) 



(7) 



with = p (f 1 2 /2 and dfl = sin 9/^6/ A <%. Ne- 
glecting all quantum mechanical interferences between 
the contributing classical trajectories the classical ioniza- 
tion probability d 3 P^/[deAdQ] = y/2e<f> Ej(^ d) ) 2 ^ 
obtained. Without angular resolution this expression has 
already been used successfully for explaining particular 
aspects of ionization by half-cycle pulses Eqs. (^) 

and (Q) constitute a new, multidimensional semiclassical 
framework within which all aspects of ionization by an 
HCP can be understood in terms of the underlying classi- 
cal dynamics of the ionized Rydberg electron. They con- 
stitute a main result of this Communication. It is worth 
noting that for all aspects which are sensitive to interfer- 
ences between different probability amplitudes, such as 
the angular distributions discussed below, the inclusion 
of the proper phases entering Eq. (ji|) is crucial. Further- 
more, it should be mentioned that there exists a wealth of 
alternative methods for solving semiclassical initial value 
problems which have been pioneered by E. Heller and 
W. H. Miller |k|. However, so far these methods have 
been applied predominantly to problems with explicitly 
time- independent Hamiltonians. 

In order to demonstrate characteristic novel phenom- 
ena in the angle- and energy-resolved ionization spectra 
and in order to show the accuracy of this semiclassical 
approach let us consider first of all the case of an HCP 
with a duration short in comparison to the classical orbit 
time of the Rydberg electron. In this case the sudden 
ionization or impulse approximation Q may be invoked. 
Thereby, the influence of the HCP on the Rydberg elec- 
tron is described as a sudden momentum change of mag- 
nitude Ap = -J^dtEit) = Ape z at t = [E(i) de- 
notes the electric field of the HCP at the position of the 
atom] . After the sudden excitation the motion of the elec- 
tron is governed by the Coulomb potential Vc(x) alone, 
i.e. A(x,i) = for t > in Eq. (|l|). Thus, immediately 
after the HCP the initial momenta of the classical trajec- 
tories appearing in Eq. (|j) are given by p^(xo) + Ap. 
The positions x^ at which the electron is imparted a 
final energy are determined by energy conservation, 
i.e. 



eo+A P 72 + pW(x 



Ap 



Af) 



(8) 



Furthermore, consistent with this approximation the 
classical actions Sj of Eq. (^) have to be replaced by 

As a general consequence of the semiclassical de- 
scription based on Eqs. (H) and (Q) scaling proper- 
ties of the underlying classical dynamics of the Ryd- 
berg electron should manifest themselves in the ion- 
ization spectra. On the classical level the trajecto- 
ries (x(t),p(t)) and (x(t),p(t)) resulting from two spa- 
tially homogeneous HCPs with field strengths E(t) and 
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E(t) = 7 4 E(7 3 i), 7 > 0, are related by the scaling 
transformation x(i) = 7~ 2 x(7 3 t), p(t) — 7p(7 3 i). In 
the impulse approximation this scaling relation has the 
far reaching consequence that the character of the clas- 
sical dynamics depends on Ap 2 and eo only through 
the ratio [Ap 2 /(— 2eo)] so that angle-resolved spectra for 
(eo,Ap, e^) and (^/ 2 e ,jAp,j 2 e^) will exhibit similar 
qualitative structure. As an example let us consider the 
case Ap 2 /(— 2eo) = 6.25 in more detail. For the sake 
of comparison with fully quantum mechanical results the 
corresponding examples are calculated for no = 50, i.e. 
Ap = 0.05 a.u., l Q = 0, and a = 0. From Eq. (g) two 
characteristic dynamical regimes can be distinguished ac- 
cording to whether is smaller or larger than the crit- 
ical energy eo + Ap 2 /2 (= 28.6 meV in this example). 
In both cases the dynamics of the ionization process are 

very different. 
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FIG. 1. (a) Full curve: initial positions of classical tra- 
jectories with no = 50, Ap = 0.05 a.u., and = 10 meV 
according to Eq. (^). The figure has to be continued into three 
dimensions by rotation around the z-axis. The numbers indi- 
cate final emission angles for specific initial positions. Angle 
counted counterclockwise from the z-axis. (b) Angular dis- 
tribution of the ionized electron ln{<i 3 Pi an / (dedf2)[a.u.]} in 
the sudden ionization approximation for the parameter val- 
ues given above: quantum mechanical result (full), uniform 
semiclassical (dashed) and primitive semiclassical (dotted). 



Classical trajectories with < e^' < eo + Ap 2 /2 start 
with negative (positive) radial momenta from the pos- 
itive (negative) z-half plane. Typical initial positions 
and final angles of trajectories which according to Eq. 
(||) yield a particular final energy are depicted in 
Fig. 0(a) for a case with e (/ ) = 10 meV. Close inspec- 
tion of Fig. |l|(a) shows that for each final emission angle 
0/ there are three different corresponding initial posi- 



tions. Therefore, one expects that the energy- and angle- 
resolved ionization probability exhibits a characteristic 
quantum mechanical interference structure which orig- 
inates from these three contributing trajectories. Fur- 
thermore, one notices that there are classical trajectories 
with 9/ = 0°(180°) whose initial positions do not lie 
on the z-axis. The existence of these trajectories and 
the axial symmetry of the dynamics around the polar- 
ization axis of the HCP give rise to the semiclassical 
glory phenomenon pT| . It manifests itself in a diver- 
gence of the classical ionization density p( d ) at 9/ = 0° 
and 9/ = 180° which originates from the coalescence 
of two axially symmeytric families of classical trajecto- 
ries. However, the resulting divergency of the primi- 
tive semiclassical result of Eq.([|) can be removed with 
the help of semiclassical uniformization methods [p|,pd|. 
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FIG. 2. (a) Same as Fig. 0(a) for e (/) = 40 meV. (b) Same 
as (a) for the pulse of finite duration described in the text. 
Shown is the region with the most significant modifications. 
Outside this region the figure would be similar to (a), (c) 
Angular distributions for the cases of Figs. [](a) (full curve) 
and (b) (dashed). 



In Fig. |l|(b) the resulting energy- and angle-resolved 
ionization spectrum is depicted. In the sudden ioniza- 
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tion approximation the result of a full quantum mechan- 
ical calculation (full curve) is compared with the corre- 
sponding semiclassical uniform result (dashed curve) and 
the primitive semiclassical result based on Eq. (J4J) (dotted 
curve). The close agreement between the semiclassical 
and exact results is apparent. The quantum interfer- 
ence between the contributing classical trajectories leads 
to characteristic oscillations of the angular distribution. 
The angular distribution is concentrated around the glory 
angles 9/ = 0°, 180° where the primitive semiclassical 
result diverges. 

The characteristic behavior of classical trajectories 
with energies 1 ^> e") > eo + Ap 2 /2 is depicted in 
Fig. §(a). In this case trajectories originating from the 
positive (negative) z-half plane start with positive (neg- 
ative) radial momenta. From Fig. ^|(a) one notices an 
extremum of 9/ on the curve of constant e^K At the 
extremum final angle 9 m the contributions of two clas- 
sical trajectories coalesce and give rise to a semiclassical 
rainbow phenomenon and a divergence of the primitive 
semiclassical approximation. Again this divergence can 
be removed by uniformization methods |j 1 1 1 . In ad- 
dition to this rainbow phenomenon glory scattering ap- 
pears at 9/ = 180°. However, the glory phenomenon at 
Qf = 0° which is present in the case of low electron ener- 
gies has now disappeared. The full curve of Fig. |(c) de- 
picts the uniform semiclassical result for d 3 Pi on /[de Adfi] 
in the sudden ionization approximation. It is again al- 
most indistinguishable from the corresponding quantum 
mechanical calculation which is not shown. For small 9 / 
the spectrum is smooth as there is only one contribut- 
ing trajectory. Around 9 m interference oscillations set 
in which reflect the appearance of two further trajectory 
classes. The three contributing trajectories persist up 
to 9/ = 180° where glory scattering is observed in the 
impulse approximation. 

Finally, let us discuss the influence of realistic pulse 
forms on these ionization spectra with the help of Eqs. 
(|J) and (0). Depending on the pulse duration and the 
spatial pulse profile various modifications of these angu- 
lar distributions are expected to appear. They will be 
studied systematically in forthcoming work ||. As an 
example, let us concentrate here on effects of finite pulse 
durations. Figure ||(b) depicts the modifications in the 
distribution of initial positions and final angles due to a 
moderately long HCP with pulse duration r = 0.01T C ; 
(T c i = 2tttiq is the classical Kepler period of the Rydberg 
electron) and pulse form E(<) = e 2 £'oexp[— (i/r) 8 ]. The 
amplitude Eo is adjusted so that |Ap| = 0.05 a.u. as in 
the previous cases. It turns out that trajectories starting 
close to the nucleus are most sensitive to the pulse shape 
of an HCP. Three major modifications in the energy-and 
angle resolved spectra can be recognized: (i) trajectories 
with initially outgoing radial momenta lead to final an- 
gles with |9/| < 37.5°, only. However, this hardly affects 
the spectrum [dashed curve in Fig. 0(c)] as the classical 



weight of this trajectory class decreases rapidly with in- 
creasing 9 f already in the impulse approximation, (ii) 
A new class of trajectories with initially incoming ra- 
dial momenta appears in the positive z-half plane. They 
have final angles |9/| < 42° and give rise to small oscil- 
lations in the spectrum at small angles, (iii) The most 
profound changes are due to trajectories in the negative 
z-half plane. Their initial positions change significantly 
so that the glory effect at 9/ = 180° disappears. Thus, 
the spectrum is most sensitive to pulse shape effects at 
large emission angles. The rainbow effect at small emis- 
sion angles is rather insensitive to the pulse shape. 

In summary, a new multidimensional semiclassical 
treatment of excitation of a Rydberg electron by a HCP 
has been presented. Based on this approach it has been 
shown that energy-resolved angular distributions of the 
ionized Rydberg electron exhibit characteristic oscilla- 
tions which can be interpreted naturally as contributions 
of interfering classical trajectories. These novel quantum 
mechanical interference features are dominated by semi- 
classical catastrophes of the glory and rainbow type. As 
all quantum mechanical phases are taken into account 
properly this semiclassical approach might also become 
a useful theoretical tool in problems concerning the re- 
construction of quantum states of Rydberg electrons with 
the help of HCPs. 
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